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Abstract 

With  the  growing  interest  in  Hall  thruster  technology,  comes  the  need  to  fully  characterize  the  plasma 
dynamics  that  determine  performance.  Of  particular  interest  is  the  existence  of  a  periodic  low  frequency 
oscillation,  commonly  referred  to  as  the  breathing  mode.  Although  there  has  been  a  significant  effort  to 
quantify  plasma  properties,  this  has  been  primarily  limited  to  time-averaged  measurements.  In  order  to 
fully  understand  the  periodic  instabilities  characteristic  of  Hall  thruster  behavior,  time  resolved  techniques 
must  be  developed.  This  study  presents  a  non-intrusive  method  of  determining  time -resolved  electron 
temperature  fluctuations  during  the  breathing  mode  cycle.  A  triggering  system  was  developed  to 
synchronize  measurements  to  thruster  discharge  current  oscillations.  Emission  in  the  NIR  is  measured 
near  the  thruster  exit  plane  using  an  optical  fiber  coupled  to  a  1.25  m  focal  length  spectrometer  and 
intensified  CCD  (ICCD)  detector.  ICCD  gating  is  controlled  using  the  output  of  the  triggering  system 
and  a  user  controlled  gate  delay.  This  allows  emission  sampling  only  from  the  desired  portion  of  the 
breathing  mode  cycle.  Emission  is  integrated  over  several  periods  to  allow  for  improved  signal  to  noise 
ratio  (SNR).  Electron  temperature  is  determined  using  the  intensity  of  Xe  I  emission  lines  in  the  context 
of  a  collisional-radiative  model.  The  resulting  measured  electron  temperature  fluctuations  during  a 
nominal  breathing  mode  cycle  are  presented  for  a  600  W  Hall  effect  thruster  at  a  single  nominal  operating 
condition. 
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Introduction 


Hall  thrusters  are  becoming  widely  used  for  spacecraft  station-keeping  due  to  their  low  thrust,  high 
specific  impulse  operation.  Despite  their  higher  system  complexity  relative  to  conventional 
monopropellant  chemical  thrusters,  the  significant  reduction  of  station-keeping  propellant  (by  as  much  as 
lOx  for  north-south  station-keeping  of  a  geosynchronous  spacecraft)  makes  this  technology  economically 
advantageous.  While  a  number  of  Hall  effect  thrusters  are  now  space  qualified,  there  is  a  strong  research 
component  actively  pursuing  the  extension  of  this  spacecraft  propulsion  technology  for  higher  powers, 
longer  life,  and  new  propellants.  These  research  efforts  have  developed  a  number  of  diagnostics  and 
models  that  presently  provide  an  adequate  understanding  of  the  time  averaged  plasma  dynamics  which 
govern  critical  processes  such  as  scaling,  sputter  induced  erosion,  and  propellant  utilization.  However, 
while  it  is  known  that  periodic  plasma  instabilities  (10’s  of  kHz  through  MHz)  strongly  affect  the  plasma 
electron  transport  within  the  discharges  of  Hall  effect  thrusters,  models  and  experiments  have  for  the  most 
part  not  measured  these  phenomena  quantitatively.  This  paper  examines  preliminary  measurements  of 
electron  temperature  using  emission  spectroscopy  with  time  resolutions  of  1  ps  over  a  periodic  15  kHz 
ionization  wave  in  the  near-field  of  a  medium  power  Hall  effect  thruster. 

The  most  often  studied  periodic  plasma  instability  of  a  Hall  effect  thruster  is  the  low  frequency  (10-50 
kHz)  fluctuation  in  discharge  current,  often  referred  to  in  the  literature  as  the  breathing  mode.  The  Hall 
thruster  breathing  mode  can  be  seen  both  experimentally  and  through  numerical  modeling1,2.  The 
oscillations  are  believed  to  be  due  to  a  cycling  of  the  neutral  and  plasma  density  in  the  exhaust  region. 
The  neutrals  fill  the  exhaust  region  and  are  quickly  ionized  due  to  the  local  high  electron  temperature. 
The  neutral  population  is  depleted  and  the  ionization  wave  moves  towards  the  anode.  The  neutral  flow 
front  moves  upstream  where  there  is  a  lower  ionization  rate.  Without  neutrals  to  ionize,  the  plasma 
density  in  the  exhaust  decreases  resulting  in  lower  electron  influx.  The  neutrals  then  refill  the  exhaust 
region,  and  the  cycle  begins  again.  This  plasma  oscillation  has  also  been  called  the  neutral  transit  time 
instability,  as  the  period  may  be  approximated  by  the  acceleration  channel  length  divided  by  the  neutral 
propellant  thermal  velocity.  However,  the  frequency  is  also  dependant  on  other  properties  such  as  the 
magnetic  field  strength,  discharge  potential,  and  plasma  density. 

Emission  spectroscopy  diagnostics  are  valuable  tools  at  characterizing  plasma  parameters,  such  as 
electron  temperature.  While  electron  temperature  has  previously  been  measured  using  various 
electrostatic  probes  3’4,  such  measurements  are  intrusive  and  have  inherent  limitations  with  regard  to  time- 
resolution.  Optical  emission  measurements  have  the  potential  to  resolve  time  scales  on  the  ns  temporal 
scale,  or  lower,  are  mainly  limited  by  SNR.  Additionally,  optical  measurements  provide  access  to  the 
thruster  near-field,  unlike  probes  which  are  limited  due  to  their  interference  with  the  plume.  Emission 
diagnostics  therefore  present  an  obvious  choice  for  the  development  of  time  resolved  plasma  diagnostics 
at  ps  or  higher  temporal  resolutions. 

The  purpose  of  this  work  is  to  demonstrate  the  capabilities  of  such  emission  diagnostics  through  time- 
resolved  measurements  of  electron  temperature  during  the  breathing  mode  cycle  in  a  medium  power  Hall 
effect  thruster  near-plume.  This  study  uses  NIR  xenon  neutral  emission  measurements  in  conjunction 
with  a  collisional-radiative  model  to  infer  electron  temperature.  Relying  on  the  highly  periodic  nature  of 
the  oscillations,  a  portion  of  the  30ps  period  is  measured  with  lps  integration  time  and  accumulated  over 
several  periods  to  achieve  adequate  SNR.  Repeating  this  process  during  the  period  of  the  breathing  mode, 
and  assuming  phase  synchronicity,  allows  us  to  measure  the  emission  temporal  behavior  over  a 
representative  period.  Therefore,  the  measurements  are  time  accurate  and  the  plasma  oscillation 
periodicity  allows  us  to  make  estimations  of  the  periodicity  of  the  electron  temperature  variation. 
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Collisional-Radiative  Model 


Theory 

A  collisional-radiative  model  (CRM)  is  used  to  relate  observed  emission  to  plasma  parameters.  A  CRM 
uses  primary  excitation  and  deexcitation  mechanisms  to  predict  emission  as  a  function  of  plasma 
conditions.  Initial  simplified  efforts  in  modeling  the  Hall  thruster  plasma  emission  were  done  using  a 
corona  model  5.  The  corona  model  simplifies  the  plasma  by  assuming  excitation  due  only  to  electronic 
collisions  and  deexcitation  due  only  to  radiative  decay.  One  difficulty  with  this  model  is  the  lack  of 
heavy  particle  collisions.  Although  electron  collisions  happen  at  a  much  higher  frequency,  it  has  been 

shown  that  heavy  particle  collisions  should  not  be  assumed  negligible  for  Hall  thruster  plasma  conditions 

6 

Another  shortcoming  of  the  corona  model  is  that  it  does  not  account  for  the  presence  of  Xe  metastable 
states.  Metastable  states  are  energy  levels  that  remain  populated  for  long  time  scales  due  to  a  lack  of 
allowed  radiative  decay  transitions,  and  therefore  collisions  with  this  metastable  state  must  be  taken  into 
account.  A  simplified  view  of  the  Xe  energy  level  diagram  with  a  number  of  radiative  transitions  of 
interest  is  shown  in  Figure  1  (a).  The  ls5  metastable  state  is  populated  from  radiative  decay  from  the  2pi 
series  but  has  no  allowable  radiative  decay  path  to  ground.  The  main  mechanism  for  leaving  this  energy 
level  is  collisional  excitation.  Therefore  this  state  remains  populated  for  longer  timescales. 


Figure  1:  (a)  Simplified  Neutral  Xenon  Energy  Level  Diagram  and  allowable  radiative  decay  pathways; 
(b)  KCD  excitation  and  deexcitation  pathways  including  metastable  (ls5)  state  and  heavy  particle 
collisions  with  ground  state 

The  KCD  collisional-radiative  model  used  in  this  study  was  developed  by  Karabadzhak  et  al7  based  on 
the  collision  cross-sections  developed  by  Chui8  and  Sommerville9.  A  schematic  of  the  KCD  model  is 
shown  in  Figure  1  (b).  The  KCD  model  builds  on  the  corona  model  by  including  heavy  particle  and 
metastable  collisions,  assumed  exclusively  with  electrons.  The  2pi  series  is  simplified  by  assuming  only 
two  possible  radiative  pathways:  radiative  decay  to  ground  through  the  ls4  state,  or  decay  to  the 
metastable  ls5  state.  Metastables  are  modeled  assuming  steady  state  equilibrium  between  metastables 
gained  through  collisional  excitation  from  ground  to  the  2pi  series  followed  by  radiative  decay  to  ls5,  and 
metastables  lost  through  collisional  excitation  from  the  metastable  state  followed  by  radiative  decay  to 
ls4. 
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The  primary  mechanisms  considered  for  the  neutral  Xe  emission  are  as  follows: 

e~  +  Xe  Xe*  +  e~  a) 

e~+Xem  -+Xe*  +  e~  b) 

Xe+  +  Xe  ^Xe*  +Xe+*  c) 

Xe2+  +Xe  ^  Xe*  +Xe2+*  d) 

The  radiation  intensity  per  unit  volume  can  then  be  approximated  using 

JxeiW  =  ^  [N0Nek*0  +  NmNekgm  +  N0Nxk^  +  N0N2k*0\  (1) 


where  Ni  are  the  species  number  densities;  k^0  is  the  excitation  rate  coefficient  for  electrons  with  neutrals 
with  emission  at  wavelength  X;  k^m  is  the  emission  excitation  rate  coefficient  for  electrons  with 
metastables;  and  k^0  and  /c^o  are  the  emission  excitation  rate  coefficient  for  Xe+  and  Xe+2  ions  with 
neutrals.  Assuming  quasineutrality,  Ne  =  N±  +  2 N2,  and  substituting  a  =  N1/Ne  ,  Equation  (1)  can  be 
rewritten  as 

Jxel  (/0  =  ~^(N0Ne)  [keO  +  a  ^10  +  ^f^2o]  (2) 

Rate  coefficients  are  determined  using 

ki  =  fiftiEdofadUidEi  (3) 

where  /,(7r,)  is  the  energy  distribution  of  species  i;  is  the  excitation  cross  sections  for  a  particle  i, 

with  energy  iq,  with  a  particle  j  to  produce  emission  at  X  ;  and  iq  is  the  velocity  of  particle  i.  Excitation 
cross  sections  for  neutral-electron  and  neutral-ion  collisions  are  determined  empirically  from  Chiu8  and 
Sommerville9.  Assuming  a  Maxwellian  electron  energy  distribution,  k^0  is  determined  as  a  function  of 
the  electron  temperature,  Te.  Rate  coefficients  k^0  and  k\o  are  determined  assuming  a  uniform  ion 
energy  distribution.  Cross  sections  for  electron-metastable  collisions  have  yet  to  be  determined 
empirically.  Therefore  the  excitation  rate  coefficient  k\m is  then  approximated  using  the  upper  2pi  level 
degeneracy  and  branching  probabilities. 


Figure  2:  Calculated  intensity  ratios  as  a  function  of  electron  temperature  using  empirical  data  for  ion 
relative  population  and  velocity  in  comparison  to  a  20%  increase  in  Xe2+  population  and  a  10%  increase 
in  ion  velocity 
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Intensity  ratios  for  the  I823/I828  and  Is34/l828  calculated  using  the  KCD  model  are  shown  in  Figure  2.  Ion 
species  populations  input  into  the  model  are  from  probe  measurements  by  Nakles  et  al10  and  Ekholm  et 
al11  and  ion  velocity  inputs  are  based  on  data  from  Nakles  et  al12.  A  parameter  variation  comparison  was 
done  using  a  20%  increase  in  the  Xe2+  population  and  a  10%  increase  in  the  ion  velocity  to  determine 
their  influence  on  the  model.  Both  models  showed  little  sensitivity  to  the  relative  species  populations  or 
ion  velocity,  except  in  the  low  electron  temperature  regions  where  ion  collisions  play  a  more  significant 
role.  The  Is34/l828  line  is  more  accurately  modeled  than  the  I823  /I828  line  ratio,  which  relies  on  the 
metastable  approximations  for  I823  line,  but  the  I834/I828  method  shows  higher  sensitivity  to  intensity  ratio 
in  the  5-10  eV  region,  making  it  more  susceptible  to  temperature  error  due  to  uncertainties  in  line 
intensity.  A  5%  uncertainty  in  the  line  ratio  in  the  5-10  eV  region  results  in  a  1  eV  uncertainty  for  the 
1-823/1828  and  a  2  eV  uncertainty  in  the  I834/I828  method.  The  lower  emission  signal  of  the  I834  line  also  makes 
this  method  more  difficult  to  accurately  determine  intensity  ratios  in  low  signal  to  noise  ratio  conditions, 
leading  to  higher  temperature  uncertainties. 


Experimental  Apparatus 
Test  Facility  and  Hall  Effect  Thruster 

The  measurements  reported  here  were  performed  in  vacuum  Chamber  1  at  the  Air  Force  Research 
Laboratory  at  Edwards  Air  Force  Base.  This  vacuum  facility  consists  of  a  2.4  m  diameter,  4.1  m  long 
cylindrical,  non-magnetic  stainless  steel  vacuum  chamber  with  two  liquid  nitrogen  baffled  (76  K),  1.2  m 
flanged  gaseous  helium  two  stage  cryogenic  (15  K)  vacuum  pumps  capable  of  pumping  speeds  of  48,500 
L/s  on  xenon.  A  cold  cathode  gauge  is  used  to  determine  a  chamber  background  pressure  of 
approximately  4  x  10"6  Torr  (corrected  for  xenon)  during  thruster  operation. 


Figure  3:  Photographs  of  BHT-600  Hall  effect  thruster.  Left:  Thruster  shown  with  optical  collection 
apparatus  mounted  to  left.  Collection  volume  taken  6  mm  downstream  of  exit  plane.  Right:  Thruster 
firing  with  Xe  propellant,  highlighting  structure  of  near  plume. 

The  Hall  effect  thruster  used  in  this  study  is  a,  600  W  BHT-600  with  a  3.2  mm  hollow  cathode 
manufactured  by  the  Busek  Company  (Natick,  MA)  This  thruster  has  been  studied  previously  using  both 
electrostatic  probes  and  various  optical  diagnostics10, 11,12,13,14,15.  Photographs  of  the  Chamber  1  apparatus 
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and  thruster  firing  are  shown  in  Figure  3.  The  BHT-600  Hall  effect  thruster  has  an  acceleration  channel 
outer  radius  of  32  mm,  an  inner  radius  of  24  mm,  and  a  10  mm  depth.  The  magnetic  field  is  produced  by 
four  outer  coils,  and  one  inner  magnetic  coil.  The  outer  and  inner  coils  are  independently  addressable  for 
optimization  of  the  field  strength.  The  nominal  conditions  and  performance  of  the  BHT-600  as  provided 
by  the  Busek  Company  are  shown  in  Table  1. 

At  nominal  conditions  the  BHT-600  produces  consistent,  periodic  40  kHz  oscillations.  A  sample  thruster 
discharge  current  measurement  is  shown  in  Figure  4.  Although  the  amplitude  of  the  signal  varies,  the 
period  of  the  fluctuations  remains  consistent  at  26.4+  0.3  ps  based  on  the  half  width  at  half  maximum 
(HWHM)  of  the  FFT  trace.  This  repeatability  is  what  allows  for  shots  to  be  accumulated  at  the  same 
position  in  the  fluctuation  over  many  periods.  Discretizing  the  current  fluctuation  into  1  ps  sections,  this 
accumulation  method  can  be  done  at  each  temporal  position  in  the  oscillation. 


Table  1:  BHT-600  Hall  thruster  at  nominal  xenon  operating  conditions  and  performance 


Parameter 

Value 

Anode  Flow  Rate 

2.45  mg/s 

Cathode  Flow  Rate 

197  (ig/s 

Anode  Potential 

300  V 

Anode  Current 

2.05  A 

Magnetic  Current  (inner  coils) 

2.0  A 

Magnetic  Current  (outer  coils) 

2.0  A 

Thrust 

42  mN 

Specific  Impulse 

1650  s 

Anode  Efficiency 

55.0% 

Figure  4:  Typical  BHT-600  breathing  mode  discharge  current  oscillations  at  nominal  conditions.  Left: 
thruster  current  oscillations  at  approximately  26  ps.  Right:  current  FFT  trace  with  oscillation  frequency 
peak  at  38  kHz. 
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Timing  and  Triggering  Circuit 

The  temporal  resolution  of  the  data  was  achieved  using  a  triggering  system  to  synchronize  the  ICCD 
gating  with  discharge  current  oscillations.  A  schematic  layout  of  the  timing  circuit  constructed  for  these 
measurements  is  shown  in  Figure  5.  The  trigger  circuit  uses  analog  components  to  process  the  discharge 
current  signal  into  a  transistor-transistor  logic  (TTL)  output,  which  is  used  to  synchronize  camera  gating 
to  the  current  oscillations.  Thruster  current  oscillations  are  measured  using  a  passive  inductive  probe, 
which  also  incidentally  acts  as  a  band  pass  filter  (3  dB  points  of  120  Hz  and  20  MHz  with  <1°  phase  shift 
at  frequencies  between  1-100  kHz).  This  low-level  voltage  signal  is  then  passed  through  a  low  pass  filter 
and  subsequently  differentiated,  so  that  the  zero  crossings  of  the  differentiated  signal  correspond  to 
minima  in  the  thruster  anode  current  signal.  Using  a  high-speed  comparator  acting  as  a  zero  crossing 
detector,  a  TTL  signal  is  produced  at  the  minimum  of  the  signal.  A  digital  delay  generator  (DDG)  is  then 
used  as  a  buffer  between  the  triggering  circuit  and  ICCD  by  not  passing  any  false  triggers  caused  by 
occasional  inflection  points  in  the  discharge  current.  The  result  is  a  TTL  signal  synchronized  with  the 
discharge  current  oscillations. 

Camera  gating  is  controlled  using  the  TTL  trigger  and  user  defined  gate  width  and  gate  delay  (down  to  ns 
timescale).  The  camera  is  triggered  using  the  TTL  output  of  the  DDG,  waits  the  specified  gate  delay,  and 
then  allows  emission  to  fall  on  the  detector  for  the  specified  gate  width.  This  timing  system  allows  the 
camera  to  gate  at  the  desired  temporal  location  over  several  periods  and  accumulate  signal  over  the 
camera  exposure  time.  The  camera  produces  a  gate  monitor  TTL  output  to  ensure  gating  is  correctly 
synchronized  to  the  desired  portion  of  the  signal.  A  1  GS/s  oscilloscope  was  used  to  monitor  the  thruster 
discharge  current,  the  DDG  trigger  output,  and  the  ICCD  gate  monitor  output  to  ensure  repeatable 
triggering  and  gating. 
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Figure  5:  Data  acquisition  timing  system  schematic 
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Optics 


The  optical  system  used  to  measure  the  emission  spectrum  is  schematically  shown  in  Figure  6.  Plume 
emission  is  collected  at  a  single  location  6  mm  downstream  of  the  thruster  exit  plane  using  a  6  mm 
reflective  fiber  optic  beam  coupler.  The  beam  coupler  focuses  a  6  mm  diameter  beam  onto  the  optical 
fiber.  The  coupler  has  greater  than  97.5%  reflectance  in  the  450  nm-2pm  range.  The  200  pm  optical 
fibers  has  a  numerical  aperture  (NA)  of  0.22  and  low  attenuation  in  the  NIR.  To  protect  the  optical 
coupler  from  redeposition  of  sputtered  material,  a  disposable  protective  1  mm  thick  quartz  window  is 
placed  over  the  optical  coupler  aperture. 

The  NIR  optimized  fiber  optic  patch  cord  from  the  optical  coupler  connects  to  a  fiber  optic  vacuum  feed 
through  the  chamber  wall  to  a  second  patch  cord  coupled  to  the  Horiba  1250M  Series  II  spectrometer 
using  a  Horiba  fiber  optic  adapter.  The  fiber  adapter  uses  100  mm  and  30  mm  lenses  to  focus  light 
emitted  from  the  fiber  onto  the  spectrometer  entrance  slit.  The  spectrometer,  with  a  1.25  m  focal  length 
(f/9),  contains  a  300  G/mm  classically  ruled  diffraction  grating,  blazed  at  500  nm,  that  produces  a 
dispersion  of  2.6  nm/mm.  For  the  measurements  presented  here,  the  spectrometer  entrance  slit  was  fixed 
at  200  pm.  While  this  limited  the  spectral  resolution,  it  increased  throughput  for  greater  signal  strength. 

The  desired  spectrum  is  detected  using  a  cooled  Andor  iStar®  Intensified  CCD  (ICCD)  with  1024  x  1024 
pixels  (13  pm  x  13  pm)  capable  of  dark  currents  as  low  as  0.065  e7pix/s.  The  camera  intensifier  system 
acts  as  a  fast  solid-state  shutter  capable  of  gating  on  1  ns  time  scales.  Gating  is  controlled  by  the  TTL 
input  and  user  defined  gate  delay  and  gate  width.  The  individual  gated  shots  are  then  accumulated  over 
the  length  of  the  desired  exposure  and  recorded  using  the  National  Instruments  Lab  View  PC  based  data 
acquisition  system. 


Figure  6:  Optical  system  used  to  measure  thruster  emission  using  a  1.25m  Horiba  spectrometer  with 
attached  Andor  ICCD.  Emission  measurements  are  triggered  using  a  custom  triggering  circuit  to 
synchronize  measurements  to  current  oscillations.  Full  breathing  mode  cycle  is  characterized  by  varying 
the  inputted  delay  between  trigger  and  ICC  gating. 
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Each  component  of  the  optical  system  introduces  a  wavelength  dependence  on  transmission.  To  correct 
for  this  detection  wavelength  dependence,  a  calibration  was  performed  using  a  NIST  traceable  200  W 
tungsten  filament  standard  of  spectral  irradiance.  The  measured  grey  body  emission  is  compared  with  the 
expected  response  to  determine  the  relative  spectral  response  (RSR)  of  the  system.  The  average  RSR  is 
then  used  to  correct  the  measured  thruster  emission. 


Results  and  Discussion 


Timing/Trigger  Circuit  Verification 

The  timing  system  demonstrated  repeatable  current  synchronized  triggering  capabilities.  A  sample 
oscilloscope  trace  is  shown  in  Figure  7.  The  2  ps  delay  between  the  minimum  of  the  current  signal  and 
the  fall  of  the  triggering  circuit  TTL  is  caused  by  the  delays  added  by  the  active  circuit  components  and 
was  found  to  be  consistent  in  all  measurements.  The  ICCD  is  triggered  after  the  initial  TTL  drop  at  time  t 
-  9  ps  and  delays  gating  until  t  -  22  ps  based  on  the  13  ps  input  gate  delay.  Additional  delay  added  by 
ICCD  electronics  was  found  to  be  on  the  order  of  nanoseconds.  After  gating  for  lps,  the  camera  waits  for 
the  next  TTL  trigger  to  start  the  next  cycle.  Signal  is  accumulated  from  each  gating  on  the  ICCD  until  the 
end  of  the  0.5  s  exposure.  Measurements  of  the  full  breathing  mode  cycle  are  achieved  by  varying  the 
gate  delay  from  0-27  ps.  Time  zero  shown  in  subsequent  figures  will  be  defined  at  the  fall  of  the 
triggering  circuit  TTL,  corresponding  to  an  input  gate  delay  of  zero.  Delays  in  gating  from  the  ICCD 
electronics  will  be  assumed  negligible.  A  +  0.2  ps  uncertainty  in  the  temporal  resolution  will  be  used 
based  on  the  HWHM  of  the  current  FFT. 


Time  (us) 


Figure  7:  Sample  oscilloscope  trace  with  ICCD  settings:  1  ps  pulse  width,  13  ps  gate  delay.  Only  80ps  of 
the  0.5  s  exposure  is  shown. 

Emission  Measurements 

A  total  of  108  emission  measurements  were  taken  at  nominal  thruster  conditions.  The  cooled  ICCD 
achieved  a  high  signal  to  noise  ratio  (SNR  >  24  dB)  for  most  cases.  Emission  measurements  near 
minimum  of  the  discharge  current  were  difficult  due  to  the  low  signal  (SNR  =1.9  dB).  Measurements 
were  taken  at  gate  widths  of  1  ps  with  gated  shots  accumulated  over  0.5  s  exposures.  The  background 
was  subtracted  from  the  raw  signal  and  the  resulting  signal  was  then  corrected  with  the  average  RSR.  A 
sample  emission  spectrum  is  shown  in  Figure  8.  A  majority  of  the  measured  lines  in  this  wavelength 
region  are  strong  neutral  lines.  The  emission  from  the  metastable  linked  823  nm  line  is  the  strongest  of 
the  three  lines.  The  metastable  state  allows  for  excited  neutrals  with  long  lifetimes,  which  can  then  be 
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easily  re-excited  producing  strong  emission  lines.  Additional  neutral  Xe  lines  seen  can  be  seen  at  820  nm 
and  827  nm,  but  their  low  signal  in  proximity  to  the  stronger  lines  makes  them  difficult  to  use  for 
diagnostic  purposes.  This  study  focuses  on  the  use  of  the  three  strongest  lines  (823,  828,  and  834  nm)  for 
determination  of  electron  temperature. 


Figure  8:  Sample  emission  spectrum  of  a  BHT-600  at  nominal  operating  conditions.  Measurement  taken 
near  peak  of  discharge  current  oscillation  with  1  ps  gated  shots  over  a  0.5  s  exposure  (19,000  cycles). 

The  peak  value  normalized  line  intensities  are  shown  in  Figure  9  with  a  representative  current  trace. 
There  is  a  strong  correlation  between  the  peak  of  the  discharge  current  and  the  peak  of  the  Xe  neutral 
emission.  Similar  results  were  seen  in  a  200  W  Hall  thruster  study  by  Liu  et  al16  using  high  speed 
imaging  to  measure  oscillations  in  visible  emission  with  discharge  current.  A  2  ps  delay  can  be  seen 
between  the  peak  of  the  normalized  intensity  and  the  peak  in  the  discharge  current.  Assuming  this  delay 
is  based  on  the  electron  transit  time  between  the  measurement  volume  and  the  anode,  an  8  km/s  average 
axial  electron  velocity  can  be  approximated.  This  is  in  approximate  agreement  with  HPHall  modeling 
which  predicts  a  5-10  km/s  electron  axial  velocity  between  the  anode  and  exit17.  A  significant  drop  in 
emission  signal  is  seen  after  23  ps,  near  the  minimum  of  the  discharge  current.  The  low  signal  in  this 
portion  of  the  cycle  increased  the  uncertainty  in  measurements,  especially  for  the  lowest  intensity  line  at 
834  nm. 


Time  (us) 

Figure  9:  Normalized  Xe  I  line  intensities  in  comparison  to  discharge  current.  Measurements  taken  6  mm 
downstream  of  thruster  exit  plane  at  nominal  operating  conditions.  Due  to  low  signal  towards  the  end  of 
the  cycle,  line  intensities  were  unable  to  be  determined  due  to  low  signal  strengths. 
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Electron  Temperature 


Electron  temperatures  are  determined  using  I823  /I828  and  I834  /I828  line  intensity  ratios  in  conjunction  with 
the  results  of  the  KCD  model.  The  resulting  electron  temperatures  are  shown  in  Figure  10.  The  I823  /I828 
method  results  in  electron  temperatures  oscillating  between  approximately  3  and  7  eV,  peaking 
approximately  3  ps  before  the  peak  in  discharge  current.  The  I834  /I828  method  results  in  larger  electron 
temperatures  fluctuations  (3-10  eV)  peaking  approximately  6  ps  before  the  peak  in  discharge  current. 
Higher  scatter  in  the  I834  /I828  method  versus  the  I823  /I828  method  is  due  to  the  I834  /I828  method’s  higher 
sensitivity  to  uncertainties  in  line  ratio.  Low  signal  at  the  end  of  the  cycle  (>  22  ps)  lead  to  higher 
uncertainties  with  the  integrated  line  intensity  for  all  three  lines,  and  therefore  resulted  in  a  higher 
uncertainty  for  both  temperature  methods. 


Figure  10:  Temporally  resolved  electron  temperature  fluctuations  measured  6  mm  downstream  of  BHT- 
600  thruster  exit  plane  at  nominal  operating  conditions.  Temperatures  are  determined  using  the  KCD 
model  a  two-line  emission  method. 

Fluctuations  in  electron  temperature  are  due  to  competing  plasma  parameters  that  are  varying  during  the 
breathing  mode  cycle.  Fluctuations  in  electron  heating  are  expected  based  on  electric  field  oscillations 
predicted  by  breathing  mode  models2.  Additionally,  changes  in  density  will  result  in  temporally  varying 
electron  collisional  cooling.  Due  to  the  system  complexities,  additional  modeling  work  is  needed  to 
determine  which  of  these  mechanisms  is  primarily  responsible  for  the  3-6  ps  delay  between  discharge 
current  and  electron  temperature. 

Also  of  note  is  the  2-3  eV  disagreement  between  the  two  methods  for  the  first  portion  of  the  breathing 
mode  cycle.  One  possible  reason  for  the  disagreement  is  the  reliance  on  the  model’s  metastable 
approximations  by  the  I834  line.  The  low  electron  density,  low  electron  temperature  segment  of  the 
breathing  mode  may  not  have  sufficient  collisions  to  maintain  the  models  metastable  equilibrium 
approximation.  Another  possible  explanation  for  this  disagreement  is  an  initial  non-Maxwellian  electron 
energy  distribution  due  to  an  insufficient  number  of  electron-electron  collisions.  Once  the  collisional 
frequency  increases  with  increasing  plasma  density,  the  electrons  return  to  a  Maxwellian  velocity 
distribution  and  are  equilibrated  by  the  increase  in  collisions. 

Experiments  to  temporally  resolve  plasma  behaviors  were  also  performed  by  Lobbia  et  al.  using  a  high¬ 
speed  dual  Langmuir  probe3.  Using  rapid  sweeping,  this  probe  is  able  to  measure  electron  density, 
electron  temperature,  and  plasma  potential  with  a  10  ps  temporal  resolution.  Electron  temperature  and 
plasma  potential  measurements  of  the  far  field  (20  cm  downstream)  of  a  600W  thruster  are  found  to  be 
nearly  in  phase  with  discharge  current.  However  the  3-6  ps  shift  seen  in  our  measurements  is 
considerably  under  the  temporal  resolution  of  the  high-speed  dual  Langmuir  probe  method. 
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Conclusions 


A  novel,  non-intrusive  method  for  time -resolved  (lps)  electron  temperature  measurement  has  been 
presented  using  emission  spectroscopy  and  the  KCD  model.  An  analog  triggering  system  was  developed 
to  synchronize  ICCD  gated  emission  measurements  with  the  desired  portion  of  the  discharge  current 
oscillation  period.  Measurements  were  taken  by  integrating  the  gated  emission  signal  over  several  cycles, 
producing  high  signal  to  noise  ratio  measurements.  The  Xe  I  NIR  emission  fluctuations  were  observed  to 
be  nearly  in  phase  with  discharge  current  oscillations.  Electron  temperature  fluctuations  (approximately 
3-10  eV)  were  found  to  be  3-6  ps  out  of  phase  with  discharge  current  oscillations.  Temperature 
disagreements  between  the  I823  /Ig28  and  Is34  /Is28  methods  in  the  beginning  of  the  cycle  may  suggest  an 
initial  non-Maxwellian  electron  energy  distribution  or  weaknesses  in  the  modeling  at  the  low  Te,  low 
electron  density  portion  of  the  cycle. 

Future  experiments  will  focus  on  achieving  spatial  and  temporal  resolution  as  well  as  expanding 
measurements  to  other  plasma  properties.  Spatial  resolution  can  be  achieved  utilizing  the  inversion 
method  presented  by  Matlock  et  al18  in  which  radial  emission  sweeps  can  be  used  to  determine  a  three- 
dimensional  estimate  of  the  electron  temperature.  Additionally,  this  optical  system  can  be  used  to 
ascertain  other  time -resolved  plasma  properties,  such  as  ion  velocity  using  Doppler  shift.  The 
development  of  these  nonintrusive  temporally  and  spatially  resolved  diagnostics  is  vital  to  fully 
understanding  the  plasma  dynamics  governing  the  Hall  thruster  breathing  mode. 
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Background:  Diagnostics 


Probe 

•  Measurements 

-  Langmuir-  Te,  density,  e-  EDF 

-  RPA-  ion  EDF 

-  Faraday-  thruster  beam  current 

•  Intrusive-  spatially  limited 

•  Temporally  limited  due  to  sweeps 

-  Lobbia1  (10  \is  resolution) 

[1]  Lobbia,  RB  and  Gallimore,  A.D. , 
Rev.  Sci.  Instrum.  81,  073503,  2010 


Langmuir  (Cold  emissive)  Probe 


Emission 

•  Measurements 

-  Line  Intensity  ratio-  Te 

-  Absolute  Intensity  -  Density 

-  Doppler  Shift  -  velocity 

•  Non-intrusive-  capable  of  near  field 
measurements 

•  Line  of  sight  averaging 

•  Measurements  on  ns  timescales 


Emission  Beam  Coupler 
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Background:  Breathing  Mode 


Constant 
neutral  flow 


o 


50 

Time  (us) 


100 


•  Seen  through  low  frequency  (10-50k  Hz)  oscillations  in 
discharge  current  (ld) 

•  Periodic  depletion  &  replenishment  of  neutrals  at  exit2  _£ 

•  Also  referred  to  as  neutral  transit  time  instability-  scales 

with  L channel  neutrals 

•  Previous  time  averaged  measurements  unable  to  ^ 

quantify  oscillations  in  plasma  properties 

[2]  Boeuf,  J.  P.;  Garrigues,  L. ,  Journal  of  Applied  Physics ,  vol.84,  no.7, 
pp. 3541-3554,  Oct  1998 
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Collisional  Radiative  Modeling  (CRM) 


*  Predicts  emission  by  modeling 
collisional  excitation  and 
allowed  radiative  decay  paths 

*  KCD3  Metastable  Modeling 

—  Treated  as  virtual  ground 
—  Assumed  in  equilibrium 

Simplified  Xe  Collisional 

Excitation  Processes 

e~  +  Xe  ->  Xe*  +  e~ 

Xe+  +  Xe  ^  Xe*  +Xe+* 

Xe2+  +  Xe  ^  Xe*  +  Xe2+* 
e  +  Xe  m  — >  Xe*  +  e 


Simplified  KCD  Diagram 


ground 


[3]  Karabadzhak  et  al.,  Journal  of  Applied  Physics ,  2006 
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CRM:  KCD  Model 
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*  Collisional  excitation  rates 

—  Ground-  experimental  cross-sections  by  Chiu  et  al.4  &  Sommerville 
et  al.5 

—  Metastable-  excitation  approximated  using  2pi  degeneracy  & 
branching  probabilities 

*  Energy  Distribution  Functions  (EDF) 

—  Ions-  assumed  uniform  velocity,  velocity  estimated  using  LIF  data 

—  Electrons-  assumed  Maxwellian  EDF,  function  of  Te 

*  Ratio  of  intensities=  f(Te),  remove  dependence  on  N0  and  Ne 

[4]Chiu  et  al,  Journal  of  Applied  Physics ,  2006.  [5]  Sommerville  et  al,  Journal  of  Prop.  &  Power,  2008. 
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BHT-600  Hall  Effect  Thruster 


•  Thruster  tested  w  Xe  at  nominal  conditions 

•  Extensive  previous  experimental  work 

-  Probe-  RPA,  Faraday,  ExB6-9 

-  Optical  measurements-LIF9’10 

[6]  Ekholm  et  al,  JPQ  2006.  [9]  Hargus  et  al,  JPQ  2008. 

[7]  Niemela  et  al,  JPQ  2006.  [10]  Nakles  et  al,  JPQ  2008. 

[8]  Nakles  et  al,  IEPQ  2009. 


BUSEK 


Dimensions 

Rj„„e,  24  mm 
Router  32  mm 
Channel  Depth  10  mm 
Nominal  Conditions 
Anode  Flow  Rate  2.45  mg/s 
Cathode  Flow  Rate  197  jug/s 
Anode  Potential  300  V 
Anode  Current  2.05  A 
Magnetic  Current  2.0  A 
Performance 
Thrust  42  mN 
Specific  Impulse  1650  s 
Anode  Efficiency  55.0% 
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BHT-600  Hall  Effect  Thruster 
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Test  Facility 


Chamber  1 ,  AFRL  Edwards  AFB 

•  Two  cryogenic  pumps 

•  1.2  m  dia.,  LN2  baffled  (70  K),  2  stage 
He  (15  K) 

•  Pump  speeds:  48,500  L/s 

•  BHT-600  operation  background 
pressure:  9.1><10'6Torr 

•  Stainless  steel  chamber  walls 

•  2.4  m  dia.,  4.1  m  length 

•  Carbon  beam  dump 

•  Flexible  1.8  mm  graphite  wall 
covering 
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System  Schematic 
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Timing  System 
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Relative  Signal 
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Relative  Signal 


Emission  Measurements 


Wavelength  (nm) 
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Normalized  Signal 


Emission  Measurements 


•  Intensity  in  phase  with  discharge  current  oscillations 

•  Similar  results  seen  in  high  speed  imaging  by  Liu  et  al.11  studying 
visible  emission  fluctuations 

•  Small  phase  shift  (~2  |xs)  seen  between  intensity  and  current 

-  Corresponding  to  a  8  km/s  electron  axial  velocity 

-  In  agreement  with  5-10  km/s  electron  axial  velocity  predicted  by  HPHall12 

[11]  Liu  et  al.,  IEEE  T.  Plasma.  Sci.,  (in  press) 

[12]  Scharfe,  Michelle  K,  Koo,  Justin  W,  personal  communication,  2011 
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Current  (A) 


Electron  Temperature 


•  Line  Ratio  Method  Disagreement  may  be  due  to  : 

-  Issues  with  metastable  approximations  in  that  region  of  the  cycle 

-  Non-Maxwellian  electron  distribution 

•  Te found  to  fluctuate  between  3-10  eV,  3-6  jjs  out  of  phase  with  ld 

-  Additional  modeling  needed  to  determine  which  heating/cooling 
mechanisms  are  causing  this  behavior 
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Future  Work 


Abel  inversion  method13 

—  Vertical  measurement  sweeps, 
Te(y),  transformed  into  Te(r) 


b(y)  =  2 


X 


r=R 


fQr^dx 


r=y 


/O')  = 


1  fR  *'00  . 

^Jr  Vy2  -  r2  y 


Full  plume  characterization,  Te(r,z,t) 

—  5  fiber  bundle  coupled  to  spectrometer, 
simultaneous  axial  measurements 

—  Use  with  Able  inversion  method  to 
fully  spatially  characterize  plume 


Compare  with  far  field  probe  measurements 


*  Demonstrate  applicability  to  other  spectroscopy  techniques 


—  Time  resolved  ion  velocity 

measurements  using  Doppler  shift 


_  V 

AA  —  H — A 


c 


[13]  Matlock,  T. ;  Larson,  C.  W. ;  Hargus,  WA;  Nakles,  MR.,  JPC,  2007 
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Summary 


*  Demonstration  of  time-resolved  (Ifis)  non-intrusive  Te  diagnostic 

9  Nominal  conditions  of  the  BHT-600  produced  3-10  eV  fluctuations  in  Te 
6  mm  downstream  of  the  thruster  exit  plane 

9  Intensity  of  NIR  Xe  I  lines  found  to  be  in  phase  with  discharge  current 

—  2  jlis  delay  associated  with  electron  transit  time  to  anode 

9  Disagreement  between  methods  for  portion  of  cycle  may  indicate: 

—  More  work  is  needed  with  the  KCD  model  metastable  approximations 
—  An  non-Maxwellian  electron  EDF 

9  Future  work  focused  on  achieving  full  characterization  of  plume, 
temporally  and  spatially  as  well  as  expanding  to  other  spectroscopy 
diagnostics 
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